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We find the admittance of a topological Josephson junction Y (ca, ifio,T) as a function of frequency 
uj, the static phase bias (fio applied to the superconducting leads, and temperature T. The dissipative 
part of Y allows for spectroscopy of the sub-gap states in the junction. The resonant frequencies 
WM,n(‘Po) for transitions involving the Majorana (AJ) doublet exhibit characteristic kinks in the 
</3o-dependence at (^o = tt- The kinks - associated with decoupled Majorana states - remain sharp 
and the corresponding spectroscopic lines are bright at any temperature, as long as the leads are su¬ 
perconducting. The developed theory may help extracting quantitative information about Majorana 
states from microwave spectroscopy. 


The interest in the condensed matter realizations of 
Majorana states is fueled by the promise of topologically- 
protected quantum computing m- While the latter 
requires the ability to braid the states, the current ex¬ 
perimental effort [SHU] focuses on indications of the 
Majorana states’ presence in various implementations of 
topologically-nontrivial superconductors. The majority 
of experiments use the dc electron transport spectroscopy 
and aim at detecting a zero-bias conductance peak asso¬ 
ciated with tunneling into a Majorana state HMg or 
the 47r-periodic phase dependence of the two “Majorana 
branches” formed by an occupied and unoccupied Majo¬ 
rana doublet [smiin]. In the former case, some addi¬ 
tional checks are necessary [5] to exclude other sources 
{e.g., Kondo effect) of the zero-bias anomaly [8l IT^ IT^ . 
Attempts to observe the unusual phase dependence rely 
on the absence of inter-branch relaxation; this is hard to 
enforce, especially over an extended time period required 
in the interference experiments [Mm ED], or at higher 
bias voltage needed for observation of multiple Shapiro 
steps |5]. 

Limitations of the techniques implemented to-date give 
an incentive to search for alternatives. Our theory eluci¬ 
dates the manifestations of the Majorana states in the 
microwave spectrum of a topological Josephson junc¬ 
tion. Spectroscopy of the Andreev (ie., sub-gap) states 
was performed recently in experiments with conventional 
metallic break junctions |2T1 E2|- The experiments did 
detect the transitions from an Andreev level to the con¬ 
tinuum of quasiparticle states [22) . and the transitions 
between the two discrete Andreev levels mi. The latter 
result in a narrow bright line in the spectrum, especially 
attractive for spectroscopy. This is why we also aim at 
a setup allowing for discrete lines in the spectrum of a 
topological Josephson junction. The junction hybridizes 
the two Majorana states to form two levels differing by 
the parity of electron number. Therefore, the particle 
number-preserving interaction with microwaves does not 
cause transitions within this doublet. That prompts us 
to consider junctions of length L > ^ allowing for higher- 


energy Andreev states, along with the Majorana doublet 
(here ^ is the coherence length in the topological super¬ 
conductor) . 

We focus on the contribution of the discrete, sub¬ 
gap states to the admittance Y{uj,(po,T) of a topolog¬ 
ical Josephson junction m|. The setup for the junc¬ 
tion is sketched in Fig. and is based on a two- 
dimensional topological insulator HZ! or a semiconduc¬ 
tor nanowire [24l [25| with strong spin-orbit (SO) inter¬ 
action in proximity with a conventional superconductor. 
The junction is controlled by the static order parameter 
phase difference (po between the leads and is probed by 
applying a small voltage V{t) induced by microwaves. It 
creates a weak time-dependent perturbation 6ip(t) of the 
phase difference, d{S(p)/dt = (2e/h)V{t), which drives 
the transitions between the sub-gap levels. The lowest 
doublet is formed by the hybridized Majorana states. 
Their crossing at is protected by conservation 

of electron number parity. (At the crossing, the number 
parity of the ground state switches between even and odd, 
see Fig. .) The crossings of the higher-energy Andreev 
levels are not protected by the parity. In a generic junc¬ 
tion, these crossings are avoided, resulting in a smooth 
dependence of the corresponding energies En{(po) on ipo, 
see Fig. 

We find the discrete lines in the microwave absorption 
spectrum originating from the transitions in the sub-gap 
energy domain. The transitions involving the states of 
the Majorana doublet, see Fig. 2a, result in a series of 
lines with a characteristic kink at po = tt in the depen¬ 
dence of the transition frequencies on ipQ, 

^TM,nipo) = E„(Tr)zfE'j^\(po-Tr\+‘^[iPo-T^)'^] ■ ( 1 ) 

Here = TdE_M / dpo are the slopes at pq = tt of the 

energies of the two Andreev levels formed by the occupied 
and unoccupied Majorana doublet at the level crossing 
point, En{'K) is the energy of the higher Andreev level, 
and the last term in Eq. ([^ represents the smooth part of 
its y>o-dependence. The intensity of these lines is quanti¬ 
fied by ReF {uj,ipQ,T), for which in the vicinity ot ip q = tt 
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FIG. 1. (Color online) Right: A typical junction setup 
hosting Majorana states. The arrowed lines mark the gap¬ 
less helical edge modes of the 2D topological insulator (2D 
TI, light grey area), proximity coupled (dashed lines) to the 
left and right superconducting leads (grey area). The junc¬ 
tion and the leads are wider than the coherence length ^ 
of superconductivity induced in the TI; ^ sets the scale for 
the bound states’ decay length into the leads. A weak time- 
dependent voltage V(t) applied between the leads gives rise to 
a small component, S(p(t) ^ 1, modulating the superconduct¬ 
ing phase difference across the junction (in the figure, the left 
lead is grounded). The resulting current response contains 
signatures of Majorana bound states. In the full setup (up¬ 
per left inset) the outer junction can be ignored since it is 
much longer then the inner o ne and gives a signal weaker by 
a factor {L/Loutev)’^, see Eq. ( |12[ ). Lower inset: the Majorana 
detection scheme is also applicable to a junction formed by a 
topological semiconductor nanowire (light grey). 


we find the following estimate: 
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Note that unlike the zero-bias anomalies in dc transport, 
the spectroscopic feature associated with the kink in the 
(^jQ-dependence of cOzfzM,n(‘Po) is not broadened by tem¬ 
perature. The brightness of lines depends weakly on tem¬ 
perature, with the scale provided by the higher-energy 
levels En- The admittance we find is a linear-response 
property of the junction; unlike the Shapiro-steps man¬ 
ifestation of the Majoranas, their effect on Y does not 
set any stringent requirement on the relaxation time of 
the system to its equilibrium state. Our main results, 
Eqs. Q and ([^, are illustrated in Fig. In the follow¬ 
ing, we outline their derivation and application to con¬ 
crete junction models. 

A microwave held induces voltage bias V{t) between 
the leads of a device built on a basis of a two-dimensional 
topological insulator (TI) or a nanowire (NW), see Fig.[^ 
The bias excites current {I{t)) between the leads, con¬ 
nected by the NW or the edge states of a TI. The ad¬ 
mittance Y{uj) of the device is dehned as a response 
function, (/(w)) = Y{uj)V{uj), at frequency oj. For the 
current operator, we may take the current through the 


FIG. 2. (Color online) The low-lying excitations spectrum 
(a) and the absorption spectrum (b) of a topological Joseph- 
son junction, (a) The dependence of the first three Andreev 
levels in the junction on the phase tpo. The lowest level has 
zero energy Em at = 71) corresponding to two decoupled 
Majorana bound states. At 930 = vr the electron number par¬ 
ity of the ground state (GS) changes. Energy Em has a dis¬ 
continuous yjo-derivative (a kink) at the GS switching point. 
The lines Ei, 2 ( 930 ) for higher levels are smooth, as the degen¬ 
eracies are lifted by disorder (^931 characterizes the avoided 
crossing). The kink in the Em{^o) function can be probed 
by microwave spectroscopy. The transitions conserving the 
electron number parity are shown by arrows 1 through 4 and 
give rise to spectroscopic lines. (Transitions 1 and 3 create a 
pair of quasiparticles above the GS.) (b) The absorption lines 
near 930 = tt, see Eq. Q. At T = 0, transitions [arrows 1 
and 3 in panel (a)] start from the GS resulting in lines 1 and 
3 (solid blue) each displaying a kink at ipo = tt. Populated 
level Em enables transitions 2 and 4 (dashed red), which also 
show a kink, c.f. Eq. §. At finite temperature this results in 
a crossing of spectral lines at 930 = tt. For the 930 dependence 
of the transitions’ oscillator strengths, see Eq. (|12[). 


junction between the TI (or NW) and the right lead, 
lit) = e - dNpi/dt. Here Npi is the number of electrons of 
the edge state in TI (or of NW) which tunneled into the 
right superconducting lead. 

We are interested in transitions between the states with 
energies below the proximity-induced gap Aq; the lat¬ 
ter inevitably is smaller than the superconducting gap 
in the leads which are the sources of proximity. That 
allows us to use the effective Hamiltonian of a proxi- 
mized TI (or NW) instead of the full Hamiltonian in 
the evaluation of dNn/dt [26]. The effective Hamilto¬ 
nian = J 111/2 in Nambu representation 

(parametrized by matrices r) takes the form 

H^°^(x) = -ihvTzUzdx - UTz + Vix)Tz -|- Mix)ax 

-k Ao(x)[t 2 , 00393 ( 0 :) — Tj, sin 93 ( 0 ;)] . (3) 


The first term here is the electron kinetic energy, fj, is the 
chemical potential, V{x) is the scalar potential (induced, 
e.g., by disorder), and Mix) is the Zeeman splitting in 
the junction; the 4x4 matrix Hamiltonian describes 
here a helical edge state in a TI [17] . (We show below that 
the effective Hamiltonian of a NW in a sufficiently large 
magnetic field takes the same form, see the paragraph 
preceding Eq. (13).) We assume the leads to be wide 
compared to the proximity-induced coherence length ^ = 
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hv/Ao [17]. That allows us to replace the leads depicted 
in Fig. by semi-infinite pads, 

Ao(a;) = Ao[0(-a:) -|- 0(a: - L)] , ip{x) = ipoQ{x - L ), 

(4) 

where L is the length of the junction. We set the order 
parameter of the left lead to be real so that (po is the 
phase difference across the junction. Using Eq. ^ to 
evaluate dNfj/dt, we find 

I={2e/h){dH^^^/dpo)- (5) 

Application of a bias V{t) to the junction results 
in a time-dependent addition Spit) to the static phase 
bias pq. According to the Josephson relation, 5p(t) = 
2eV{t)/h. For a monochromatic bias voltage V{t) = 
V{ui) cos ujt, we can treat Spit) perturbatively as long 
as \eV{uj)/huj\ 1. (Hereinafter, we set h = 1.) In this 
weak bias limit, the Hamiltonian H can be split into a 
time-independent part and a time-dependent per¬ 

turbation 


Idip^). (6) 

To evaluate the admittance Yipj), we apply the standard 
linear response theory to the problem set by Eqs. (i-®. 
The result may be expressed in terms of the spectrum 
of and the matrix elements 

J dx^*^{x)[d^„n^°\x)]<i>n{x) (7) 

of the operator defining the perturbation, see Eq. ®. 
Here $„(a;) are the eigenfunctions of and energies 
En are measured from the Eermi level. As the result of 
using the Nambu spinor notation, the energy spectrum is 
particle-hole symmetric (PHS), meaning that the eigen¬ 
values come in pairs {En, —En)- We will label by —n the 
state with energy —En for all n > 0; the lowest (n = 0) 
doublet is {Em, —Em)- We concentrate on the absorp¬ 
tion lines which originate from transitions involving a 
state of this doublet. Using the Kubo formula [28] and 
PHS, we obtain for the corresponding part of the admit¬ 
tance: 

R«rM = — E 

^ E„>Em 

^ - ^-M,n) + S{huJ - UJ+M,n)] - 

( 8 ) 

Here uJipM.n = En T Em, see Eq. ®. Note that there 
are no terms with En = Em since n^._M = 0 by Pauli 
exclusion principle. The eigenvalues En with n ^ 0 are 
in general not degenerate, contrary to the case of a con¬ 
ventional time-reversal symmetric (TRS) S-N-S junction, 
where the states are doubly degenerate (Kramers dou¬ 
blets). The “degeneracy” and zero value of the energy 


, 1 En 1 ^ E-n 

tanh-tanh- 

2T 2T 


Em sX pq = e are protected by symmetries, and will be 
discussed below. 

The junction sub-gap excitation spectrum is found 
from Eqs. ® and Q by using the standard scattering 
matrix method . We find that there is always a state 
with vanishing energy Em sX pq = e [26] . The Em = 0 
state is protected by fermion number parity conserva¬ 
tion. (We assume the junction is well separated from 
other junctions or interfaces with bound states, allowing 
us to ignore hybridization of those states with Em, see 
Fig. a and footnote m-) Near po = e the dispersion 
Em{‘Po) is linear. For a reflectionless junction [M = 0 in 
Eq. ®] we find 

1 V 

EM{‘fo) = ^j^^\^o-E\- ( 9 ) 

Breaking of TRS in the junction, M ^ 0, results in 
backscattering. This does not lead to qualitative changes 
to Eq. § but merely modifies the prefactor in it |30j . 

The situation is different for the higher Andreev lev¬ 
els. In a reflectionless junction they have degeneracies 
at ipo = e: in our notation, for odd n the levels n and 
n-l-1 are degenerate, En\E) = £'i+i(7r) = (n-|-l)7™/2L. 
(We take now L ^ Backscattering lifts these de¬ 
generacies and thus leads to qualitative changes in the 
spectrum dX = e (see Fig. [^). The resulting avoided 
crossing makes the :^o"dependence of the levels’ energies 

smooth m 

En{po) = eI^\e) + (-l)”^\/(‘^o-7r)2 + (5(^„)2. 

( 10 ) 

(We ignore here corrections to the prefactor v/2L due 
to weak backscattering.) The smoothness in ( |Io| ) is 
characterized by the width of the avoided-crossing re¬ 
gion 5pn = \r+ We denote by r^/_(E) 

the reflection amplitudes for electrons/holes entering the 
junction from the left (right) at energy E- In a sim¬ 
ple model with U(a:) = 0, M{x) = MQ{x)Q{L — x), 
and M ^ v/L, the junction is symmetric and r±{E) = 
M(l —± p,). In this model, for small 
chemical potential, p <C v/L, the width of the avoided 
crossing is [32] 


5pn = 


\M\ Ia^I 
^^/2FF;i°)(7r) 


( 11 ) 


The sub-gap spectrum determines the resonance fre¬ 
quencies u}^M,n{'Lo)' combining Eqs. ® and (10) yields 
Eq. ([^. The brightness of the spectral lines is set by 
the matrix elements in Eq. ®. They are calculated from 
Eq. ® where, according to Eq. Q, the integration over 
X is restricted to the right superconducting region x > L 
where d,pg'H^^\x) is non-zero. There, the sub-gap wave 
functions decay exponentially, $„(a;) = e~^^~^'>/^^n{L)- 
We assume here En ^ Aq so that the decay length is 
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approximated bv £■ In the limit of weak reflection, we 
find [26] 


lO/lb |2_ ^ 
\^±M\n\ ~g]^2 


i±(-iy 


, lyjQ - 7r| ± 2Re6^*j^S(pne^'^^ 


( 12 ) 

Importantly, 1/^^; likewise, the contribu¬ 

tion in the admittance from the long outer junction in 
Fig. 0 is proportional to and can be ignored 

since the corresponding length Fouter ^ L. In Eq. 

■dn = and 5<^m = ii[d\dEr%]E=Q 

is the small correction (due to backscattering) to the wave 
function of level Em\ here d± is the transmission am¬ 
plitude. In the simple model used to derive Eq. 0 , 

Since Eq. (12) was derived assuming L ^ we can 
express its prefactor as (£'(y^)^/2 by using Eq. § in that 
limit; the replacement is valid if the levels resolved in the 
transitions are far below Aq. Furthermore, close to tpo = 
TT we can neglect the second term in the square brackets 
in ([Ill) as long as |(^o — 7 r| ^ since |<54)^| «C 1 . 

Using the approximate matrix elements in Eq. (jsj) leads 
to Eq. ([ 2 ]). 

Approximating the matrix elements in Eq. (§ by a 
constant (if (^)^/2 does not capture the brightness vari¬ 
ations of the absorption lines with n and pq. At (/jq = 
transitions to levels En with odd n (lines 3 and 4 in 
Fig.Uja) are brighter than transitions to levels E^+i (lines 
1 and 2 ), due to <5$^ 7 ^ 0. As a function of tpo, the 
brightness is non-analytic aX pq = tt which is a maximum 
(minimum) for lines 1 and 4 (2 and 3) in Fig. Far 
from the avoided crossing, |(/?o — 7r| ^ dpn, the spectral 
lines 1 and 4 become dim. In general, the lines corre¬ 
sponding to frequencies ui-M,n and co+M.n+i with n odd 
become dim away from po = tt. This is because of an ap¬ 
proximate selection rule for the matrix elements 
we see from Eq. ( jlij ) that and for 

odd n are smaller by factor oc 'R,e5^*j^5ipnE^''/\po — t^\ 
compared to those with even n. |33j This analysis allows 
us to extrapolate our theory to stronger backscattering, 
dpn ~ 1. Due to the large width of the avoided crossing, 
the above approximate selection rule becomes inapplica¬ 
ble. On the other hand, at pQ = e the alternation of 
the lines’ intensities with n becomes more pronounced. 
The main feature, the kink in the transition frequencies 
at ipo = TT, persists. 

We derived Eqs. 0 and 0 , for a TI junction, but the 
same low-energy model, Eqs. 0 -([To|), is applicable to a 
NW-based setup. For illustration, we concentrate here 
on the limit of large Zeeman energy, B ^ ma^] the SO 
energy scale here is determined by the electron effective 
mass m and SO velocity a. At low-energies E B, 
linearization |34| of the spectrum near the Fermi points 
k « ±kz = ±\/2mB leads to an effective Hamiltonian of 
the form 0 with v = kzlm^ and Ao = A2ma/fcz where 


A is the induced s-wave gap in the nanowire. (We set 
fj, = 0 for simplicity.) The structure of spectrum of a 
long NW junction is therefore identical to that of the TI 
junction described above - only the microscopic forms of 
the phenomenological parameters in 0 are different. 

In the NW junction even TRS scalar disorder may lift 
the degeneracies of high Andreev levels. We illustrate 
this by considering a short-range impurity uo6{x — L) in 
the microscopic Hamiltonian. (For definiteness, we take 
the impurity to be at the junction interface, as was done 
in Ref. I24l j The low-energy projection of the impurity 
Hamiltonian yields respective forward and backscattering 
terms [22] V{x) « uoS{x — L) and M{x) « uq5{x — L) in 
Eq. 0. (The prefactors here are given to lowest order 
in /B.) The width of the avoided crossing. 


5pn 


|uo| g Ei°^(7r) 

V L v/L 


(13) 


is finite due to M yf 0; here, as before, for n odd 


Ei°^(7r) = ifl'((^(7r) = (n -|- l)7ru/2(L -|- ^). The correction 


p(0) 


(7r) = F;),+iv 


to the lowest-level wave function is = iuo/v){l — 

iuo/v) and the phase appearing in Eq. (12) is = 
( mo /| uo |)(1 — iuo/v). We see that the avoided crossings 
and the transition matrix elements in a NW junction can 
be quantified in the same way as in a TI junction. 

Using Eq. (9), we may compare the strength of ab¬ 
sorption, Eq. (2), with the corresponding strength |35| of 
transition within the pair of Andreev levels in a conven¬ 
tional short S-N-S junction 1211112]. At equal frequencies 
w, the transition in the “Majorana” junction is stronger 
than the one in a conventional junction if the transmis¬ 
sion coefficient of the latter is < 0.4. 

In summary, we have shown how Majorana bound 
states manifest in the finite-frequency admittance of a 
topological Josephson junction with multiple Andreev 
levels. Our main finding is a kink in the (/?o"dependence 
of the resonant absorption frequency, and can be ob¬ 
served in the dissipative (real) part of the admittance 
y(w), see Eq. 0 and Fig. 0. Alternatively, one may 
employ the reflected microwaves’ phase shift, argy(a;), 
obtained from Eq. 0 and its Kramers-Kronig partner. 
The frequency-dependent phase shift jumps by tt across 
a resonance, and the position of this jump as a func¬ 
tion of ipo shows a kink. The kink is a consequence of 
the ground state parity switching in the junction, or de¬ 
coupling of Majorana states. The admittance provides a 
novel, linear-response signature of this decoupling. 

We thank Michel Devoret and Liang Fu for discussions, 
and Richard Brierley and Hendrik Meier for valuable 
comments on the manuscript. This work was supported 
by NSF DMR Grant 1206612, ONR Grant Q00704, ARO 
Grant W911NF-09-1-0514, DFG through SFB 767, EU 
FP7 Marie Curie Zukunftskolleg Incoming Fellowship 
Programme (Grant 291784) and the Senior Fellowship 
of the Zukunftskolleg, Konstanz. 












5 


[1] A. Kitaev, Annals of Physics 303, 2 (2003) 

[2] C. Nayak, S. H. Simon, A. Stern, M. Freedman, and 
S. Das Sarma, Reviews of Modern Physics 80, 1083 
(2008), arXiv:0707.1889 [cond-mat.str-el] 

[3] J. Alicea, Reports on Progress in Physics 75, 076501 
(2012), arXiv:1202.1293 [cond-mat.supr-con] 

[4] C. Beenakker, Annual Review of Con¬ 
densed Matter Physics 4, 113 (2013) 

http: //dx.doi.org/10.1146/annurev-conmatphys-030212- 
184337 

[5] V. Mourik, K. Zuo, S. M. Frolov, S. R. Plissard, E. P. 
A. M. Bakkers, and L. P. Kouwenhoven, Science 336, 
1003 (2012) 

[ 6 ] L. P. Rokhinson, X. Liu, and J. K. Furdyna, Nature 
Physics 8 , 795 (2012) 

[7] M. T. Deng, C. L. Yu, G. Y. Huang, M. Larsson, 
P. Caroff, and H. Q. Xu, Nano Letters 12, 6414 (2012) 
pMID: 23181691, http://dx.doi.org/10.1021/nl303758w 

[ 8 ] A. Das, Y. Ronen, Y. Most, Y. Oreg, M. Heiblum, 
and H. Shtrikman, Nature Physics 8 , 887 (2012) 
arXiv:1205.7073 [cond-mat.mes-hall] 

[9] G. Koren, T. Kirzhner, E. Lahoud, K. B. Chashka, 
and A. Kanigel, Phys. Rev. B 84, 224521 (2011) 
arXiv:1111.3445 [cond-mat.supr-con] 

[10] C. Kurter, A. D. K. Finck, Y. S. Hor, and D. J. Van 
Harlingen, ArXiv e-prints (2013), arXiv:1307.7764 [cond- 
mat.mes-hall] 

[11] S. Hart, H. Ren, T. Wagner, P. Leubner, M. Miihlbauer, 
C. Briine, H. Buhmann, L. W. Molenkamp, and A. Ya- 
coby. Nature Physics 10 , 638 (2014), arXiv:1312.2559 
[cond-mat.mes-hall] 

[12] V. S. Pribiag, A. J. A. Beukman, F. Qu, M. C. Cassidy, 
C. Charpentier, W. Wegscheider, and L. P. Kouwen¬ 
hoven, ArXiv e-prints (2014), arXiv: 1408.1701 [cond- 
mat.mes-hall] 

[13] K. T. Law, P. A. Lee, and T. K. Ng, Phys. Rev. Lett. 
103, 237001 (2009) 

[14] K. Flensberg, Phys. Rev. B 82, 180516 (2010) 

[15] J. D. San, S. Tewari, R. M. Lutchyn, T. D. Stanescu, 
and S. Das Sarma, Phys. Rev. B 82, 214509 (2010) 
arXiv:1006.2829 [cond-mat.supr-con] 

[16] H.-J. Kwon, V. M. Yakovenko, and K. Sengupta, Low 
Temperature Physics 30, 613 (2004) 

[17] L. Fu and C. L. Kane, Phys. Rev. B 79, 161408 (2009) 

[18] H. O. H. Churchill, V. Fatemi, K. Grove-Rasmussen, 
M. T. Deng, P. Caroff, H. Q. Xu, and C. M. Marcus, 
Phys. Rev. B 87, 241401 (2013) 

[19] A. D. K. Finck, D. J. Van Harlingen, P. K. Mohseni, 
K. Jung, and X. Li, Physical Review Letters 110, 126406 
(2013), arXiv:1212.1101 [cond-mat.mes-hall] 

[20] D. M. Badiane, L. 1. Glazman, M. Houzet, and J. S. 
Meyer, Comptes Rendus Physique 14, 840 (2013) 

[21] L. Bretheau, (/. Girit, H. Pothier, D. Esteve, and 

G. Urbina, Nature 499, 312 (2013) 

[22] L. Bretheau, (/. O. Girit, C. Urbina, D. Esteve, and 

H. Pothier, Phys. Rev. X 3, 041034 (2013) 

[23] We consider a system where the phase ipo is classical. 
Our method does not require its quantum fluctuations, 
unlike, for example, the proposal in Ref. m- 

[24] R. M. Lutchyn, J. D. San, and S. Das Sarma, Phys. Rev. 
Lett. 105, 077001 (2010) 


[25] Y. Oreg, G. Refael, and E. von Oppen, Phys. Rev. Lett. 
105, 177002 (2010) 

[26] See Supplemental Material for details. The supplement 
includes references to iJTHini 

[27] In the case of a NW junction, we assume that the ends of 
the wire are much farther than ^ away from the junction. 
In the TI setup there are two junctions, see Fig. The 
outer one can be ignored if it is much longer than the 
inner one, see remark below Eq. ( |12[ ). 

[28] L. Landau and E. Lifshitz, Statistical Physics, 3rd ed., 
Vol. 5 (Butterworth-Heinemann, 1980). 

[29] C. W. J. Beenakker, Phys. Rev. Lett. 67, 3836 (1991) 

[30] The slope is dEM/d(fio = -\- 0{\r+\^)), see 

Ref.IM 


[31] Equations 0 and ( |10[ | are valid in the full range of (po ex¬ 
cept for a narrow interval around <po = 0 where a similar 
avoided crossing happens. 

[32] Lifting of degeneracies in the model for which Eq. 0 
was derived requires, in addition to M yf 0, a non-zero 
chemical potential p, due to special symmetries [26] of 0 
at the Dirac point. 

[33] This selection rule is exact in a reflectionless junction, 
5 v5n = 0, where the operator commutes with and 

hw. 


[34] J. Klinovaja and D. Loss, Phys. Rev. B 86, 085408 

( 2012 ) 

[35] F. Kos, S. E. Nigg, and L. 1. Glazman, Phys. Rev. B 87, 
174521 (2013) 

[36] E. Ginossar and E. Grosfeld, Nature Communications 5, 
4772 (2014), arXiv:1307.1159 [cond-mat.mes-hall] 

[37] C. Wang, Y. Y. Gao, 1. M. Pop, U. Vool, C. Axline, 
T. Brecht, R. W. Heeres, L. Frunzio, M. H. Devoret, 
G. Catelani, L. I. Glazman, and R. J. Schoelkopf, Nature 
Communications 5, 5836 (2014) 

[38] A. C. Potter and P. A. Lee, Phys. Rev. B 83, 184520 

( 2011 ) 

[39] C. W. J. Beenakker, D. 1. Pikulin, T. Hyart, H. Schome- 
rus, and J. P. Dahlhaus, Phys. Rev. Lett. 110, 017003 
(2013) 

[40] S.-f. Zhang, W. Zhu, and Q.-f. Sun, Journal of Physics: 
Condensed Matter 25, 295301 (2013) 









Supplemental Material 


Jukka 1. Vayrynen,^ Gianluca Rastelli,^’^ Wolfgang Belzig,^ and Leonid I. Glazman^ 

^Department of Physics, Yale University, New Haven, CT 06520, USA 
^ZukunftskoUeg, Universitdt Konstanz, D-78457 Konstanz, Germany 
^Faehbereich Physik, Universitdt Konstanz, D-78457 Konstanz, Germany 


In this supplemental material we show in detail the derivations that were omitted in the main text. 


SMI. DERIVATION OF THE CURRENT OPERATOR 


To model the S-TI-S (or S-NW-S) junction, we assume that the left (x < 0) and right (x > L) halves of the 
one-dimensional edge are tunnel-coupled (with tunneling Hamiltonians and Hj^t) to the respective leads. The 
junction comprises of the segment 0 < x < L, where L is the distance between the leads, see Fig. 1 of the main text. 
The operator of tunneling current into the right lead is 

/ = (SI) 

Since Hut conserves the total number of electrons, it commutes with Njn, + Nfi, where Nup and Nn are the 
respective electron number operators for the right lead and the part of the edge in tunnel-contact with the lead [SI]. 
The conservation law allows us to write Eq. (SI) as / = —e{i/H)[HRT, Hr], Since we are interested in the physics 
of the junction bound states, it is beneficial to move to a low-energy description that includes only the edge degrees 
of freedom. This is possible because the superconducting leads are gapped; the description is valid at energies less 
than the pairing gaps of the leads, El <C A. The current operator / projected to the low-energy degrees of freedom 
is easiest found by writing it first in a form that does not explicitly contain the lead operators. Denoting H' the full 
(bare) Hamiltonian of the edge, and H'j^ its kinetic energy term, we can write [Hrt, A^i?] = [{H' — H'j^), (since 
the only terms in H' that do not commute with Nr are H'j^ and Hrt). This yields 

I = -e{i/n)[{H' -H'i,),Nr]. (S2) 

Now we are ready to project / to the low-energy subspace. At low energies E ^ A, single electrons from the edge 
cannot tunnel into the leads. Pairs of electrons, however, can tunnel, and this gives rise to an effective pairing between 
them. When moving to the low-energy description, the aforementioned tunneling term Hrt is replaced by a particle 
non-conserving pairing Hamiltonian H/^ r, given by 


-y i'L-\-a 

Ha,r = -rAo dxA/^Tx cos (p - Ty sin , 

L 


(S3) 


where a is the length of the right proximitized region (see Fig. 1 of the main text), and (p is the relative phase 
between the leads (we set the phase of the left superconductor to zero). The magnitude of the proximity induced 
gap at energies E <C A is Aq = ^jA/{A 4- ^ 7 ) [S 2 ], where 7 is the tunneling rate between the edge and the lead in 
absence of A. In Eq. (S3) we introduced the spinor notation = (f/jj, ipi, where creates a right 

(left) moving electron on the edge. The Pauli matrices Tx,y,z act in the particle-hole space. In the spinor notation 
Nr = dxA>^TzA! /2. Projecting the current operator I to low energies, we can replace {H' — H'j^) —> in Eq. 

(S2). Gommuting Nr with Ha,r, we find that the projected current operator can be written as 


. ^ ^dH 


(S4) 


since Ha.r is the only i^-dependent term in H. 

To irradiate the junction with microwaves, one applies a weak time-dependent voltage V(t) to, say, the right 
lead, where we also measure the current. (The specific relation between V{t) and the electric field amplitude of the 
microwaves depends on a concrete device geometry [S3].) The effect of the resulting bias voltage is to make the low 
energy Hamiltonian H time-dependent: the phase difference p in Eq. (S3) has to be replaced by p{t) = 4- Sp{t), 

where po is the phase difference in absence of bias and Sp{t) satisfies the Josephson relation Sp{t) = 2eV{t)/h. The 
replacement can be derived by doing a time-dependent gauge transformation that removes the bias from the leads. 
The transformation gives rise to a time-dependent phase in the, say, right tunneling amplitudes and, at low energies, 
in the order parameter in Eq. (S3). Developing (S3) to first order in Sp{t) leads to Eq. ( 6 ) of the main text. 
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SM2. THE FULL EXPRESSION FOR THE ADMITTANCE 


The admittance Y{uj) is found by using the standard Kubo formula with H^) [Eq. (6) of the main text] as the 
perturbation. We find (a;’*' = w + Of), 




+ 


uj~^Lj 

8ie^ 




E 


E„>Em>0 


(w+)2 - (E„ - EmY 


IfiE^) - /(£;„)] 


(w+)2-(E„+E^)2 


[1 - fiEm) - fiEn)] ■ (S5) 


Here f{E) = l/(e^/^ + 1) is the Fermi function. The transitions in (S5) conserve the fermion (and quasiparticle) 
number parity: the first term in the brackets arises from transitions between two different quasiparticle states, while 
the second term accounts for creation/annihilation of a pair of quasiparticles. There are no terms that create an 
odd number of quasiparticles. Also, there are no terms with Em = En since none of the levels are degenerate in a 
generic disordered junction and 'H^m-n = 0 by Pauli exclusion principle. The first term in Eq. (S5) represents the 
conventional inductive response of a Josephson junction. The inverse inductance 1/Lj = 2ed^glj is a derivative of the 
Josephson current, Ij{(po) = {I)o, with the current operator defined in Eq. (5) of the main text, and ■ ■)o denoting 
the averaging over an equilibrium density matrix with the Hamiltonian (3) of the main text. 

Equation (8) of the main text is obtained from Eq. (S5) by taking the real part and accounting only for the positive 
frequency transitions involving the lowest level, Em = Ej^ ■ 


SMS. THE EQUATION FOR THE BOUND STATE ENERGY SPECTRUM 

The bound state energy spectrum is calculated by using continuity of the eigenstate wave function $£;(x) at the 
two interfaces, x = 0, L. This leads to the condition (see Sec. SMS A and Ref. S4) 

det[l - Sn{E)Sa{E)] = 0, (S6) 

which determines the energy spectrum. Here Sm and Sa are respectively the normal and Andreev scattering matrices 
(defined in Sec. SMS A). The former describes the reflection in the junction (0 < x < L), while the latter accounts 
for the Andreev reflection at the two interfaces at x = 0, L. We assume hereafter that there is backscattering only in 
the junction, i.e., the Zeeman field is vanishing outside it, M(x) = 0 for x < 0 and x > L [S5]. In particle-hole space 
the two matrices are 



Sa{E) = 


0 

0 ) 


(S7) 


where e = /S.q/[E + inv), and kv = yj Aq — E"^ is real for bound states. The components in spin-space are 



r±(E) d'^{E) 
d±{E) r'^{E) 


, ta = 


1 

0 



(S8) 


The scattering amplitudes r^{E), d^^{E) in terms of the junction parameters are given in Sec. SM4. Hereinafter, 
we omit their energy argument in the notation, unless there is a risk of confusion. Inserting Eqs. (S7) and (S8) into 
Eq. (S6), and using the unitarity of S± yields the equation for the spectrum. 


^ 2 ^a(E) - r+r\){d-d'_ - 

-d+d'_e-^^» - d'^d_E^° - r+r_ - r'^r'_ = 0 . 


(S9) 


The equation (S9) is difficult to solve in full generality, and its tractable limit will be considered in Sec. SM6. However, 
already at this stage we note that this equation has a solution Ea 4 {ipq) that vanishes at (po = tt. This can be seen 
by using in Eq. (S9) the following two properties of the scattering matrices S± [which follow from two symmetries 
of the Hamiltonian (3) of the main text]. First, from particle-hole symmetry, one has S-{E) = —azS^{—E)az [see 
Eq. (S42); * denotes complex conjugation]. Second, from the invariance of junction Hamiltonian under the combined 
reversal of time and magnetic field [M(x) —)■ —M(x)], one finds d± = d'^ [see Eq. (S44)]. These two properties of 
the scattering matrix are derived in Sec. SM4B. Using the two properties as well as the unitarity of the scattering 
matrix, Sl.S± = I, it is straightforward to check that Eq. (S9) has a solution Ea 4 = 0 at (/?o = tt. 
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By expanding Eq. (S9) in (/?□ — tt ^ 1, we can evaluate the slope of the lowest Andreev level. We find 

(see Sec. SM3B), 

(SIO) 

^\dEd+ + id+Ao Y - \dE{d+r*_)Y 

where the scattering amplitudes on the right-hand-side are evaluated at A = 0 (in which case |(i+| = |d-|). Let us 
consider two simple limits of Eq. (SIO). First, when Aq is much smaller than the energy scale at which the scattering 
amplitudes vary, we recover ^^pgEJ^(^T) = ±|(i+|Ao/2 which was already found in Ref. [S5] for a junction of length 
L ^ with M{x) = M<d{x)Q{L — x). The second example is that of a weak magnetic field M{x). We can then 
use the Born approximation (BA) for the scattering matrix. The reflection amplitudes are small, |r£^| ^ 1, and the 
second term in the denominator of Eq. (SIO) can be neglected. In the model of Eq. (3) in the main text, the BA is 
valid when dxM{x)/v <C 1, and amounts to setting d+ = 6*“^+ with dEd+ = L/v (see Sec. SM4). Equation (SIO) 
then becomes 

9^„i?Ar(7r) = ±i^+0(|r+n, (Sll) 

leading to Eq. (9) of the main text (see also the footnote 30 there). The terms (!l(|r+p) are given in Eq. (S88) in 
Sec. SM6E. Finally, we wish to emphasize that Eqs. (SIO) and (Sll) are true for a generic disordered junction with 
any scalar disorder potential V{x) that does not exceed the TI bulk band gap. The problem is tractable, because 
of the topological protection of the helical states against scalar disorder - the potential V{x) does not backscatter 
electrons and the junction is reflectionless in absence of M{x). 


A. Derivation of Eq. (S6) and the wave function boundary values 


We write the wave function as <!> = (u^, ui, vi, v^)’^, suppressing the energy labels in this section. The scattering 
matrices for the two superconducting regions (x < 0 and x > L) are defined as 


M^(-OO) 

u^(0) 

x^(-oo) 

Vt(0) 


M^( —oo) 

Wt(0) 


Ui{L) \ 

u^(oo) _ 


x^(-oo) ’ v^{L) v^{L) 


u^{L) 

M;(oo) 


w^(oo) 


where the incoming states are on the right. (We shall shortly require $(±00) = 0 since we are considering bound 
states.) The Andreev scattering matrix for the left/right “semi-infinite” superconducting region is 


Sa,lir = 


qI’Vl/r 


q—i-Vl/r 


^—ia{E) _ ^0 


E + ivn’ 


= -JaI-E\ (S13) 


and is derived in Sec. SM4A. For bound states $(±00) = 0, and Eq. (S12) can be replaced by the single equation 
/ wt(0) \ / M;(0) \ 


wt(0) \ / M;(0) \ 

^t(o) " ^ n(o) ’ 0; ’ 

Vi(L) J V v^{L) J 


( 0 

ta \ 



0 \ 


0 

II 

[0 

g*V0 J 


where we set (pL = 0, tpR = ipo. 

In the junction, 0 < x < L, we have the normal state scattering relation 


wt(0) \ / w^(0) \ 

UtW _ q Ul{L) q _ f S+ 0 \ 

n(0) \ v,{0) h 0 S.)’ 

v^{L) J V vAL) J 


r± d± 
d± r± 


where and are the reflection and transmission amplitudes for electrons/holes entering the junction from 

the left (right). The eight equations in (S14) and (S15) allow us to express all the components of <I)(0) and $(L) in 
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terms of a single component, say v^{L), which can be taken to be the wave function normalization constant. One of 
the eight equations remains unused - this is Eq. (S6) which determines the energy spectrum. Explicitly, it reads. 


- d+{d-d'_ - - d-{d+d'+ - 

— (d+r_ + {d-r+ + = 0 . 


(S16) 


Using unitarity of the scattering matrices, S'±S'± = 1, in (S16) leads to Eq. (S9). 
The wave function components are 


and 


u^(0) 


n(o) 


Wi(0) 


ut(0) 


rV(d-d'_-r_/_)e—(^)+r-e-(^) 

_ d+(d_d'_ _ r_r^_)e“*“(®)e“®‘^o ^ ’ 

r'+{d-d'_ - r_r'_)e-“(^) + 
c?_e*“(^) _ d+((i_d'_ _ ^ ’ 

rL(d+d; - r+r'+)e-“(^) + 
d_r+ + d+r'_e~'^^° 
r^d+dV - r+rV)e—(^) + r+e-(^) 

d_r+ + d+r'_e~'-^° ^ 


m^(L) 




_ d+r_ + d_r+e*‘^« 

d_e®“(^) — d+(d_d'_ — 

_ d+r_ + d-r'^e^‘^° 

d_e®“(^) — d+(d_d'_ — 




(S17a) 

(S17b) 

(S17c) 

(S17d) 

(S17e) 

(S17f) 

(S17g) 


The spatial dependence of <I)(a:) can be found by acting on, say, <f’(L) with the transfer matrix, see Eq. (S31) in 
Sec. SM4. 


B. Derivation of Eq. (SIO) from Eq. (S9) 

It is useful to divide Eq. (S9) into two parts, 

/i(E) + /2(E,^o)=0, 

where /i does not contain explicit (/SQ-dependence, 

/i(E) = + (d^ - r+r'+){dl - r_r'_)e-2“(^^) 

f2iE,(po) = -2d+d_ cos(,?o • 


r+r_ — r'^r'_ 


(S18) 


(S19) 


We have used here the property d = d^, see Eq. (S44). Using d/dipo = d^o^^E + d^o’ the hrst derivative of Eq. (S18) 
yields {E' = d^gE) 


E'dEfiiE) + E'dEf2{E,^o) + d^gf2{E,^o) = 0 , 


(S20a) 


Note that at = tt we have d,^„f 2 = 0 and thus Oe/i = —dEf 2 from Eq. (S20a). Taking the second derivative of 
Eq. (S18) yields 

E'2 (a|/i(U) + d%f 2 {E, ^o))+E" {dEfi{E) + dEf 2 iE, ^o))+E'd^gdEfi{E)+2E'dEd^gf2{E, (^o)+9^„/2(E, ^o) = 0 . 

(S20b) 

Evaluating this at (^o = '^ and using Eq. (S20a) leads to 


{d%fi{E) + dlf2{E, <^o)) + dlgf2{E, (^o) = 0 . 


(S21) 
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From Eq. (S18) we have /i = —/2 = 5^„/25 and thus 


{d^Ef = 


dl{h{E) + ME,7T)) 


Next we will evaluate the derivatives in the denominator. From particle-hole symmetry, Eq. (S42), we have the 
following properties valid at zero energy: = (—and = (—for the nth derivative, 

n = 0 meaning no derivative. Using Eq. (S19) and the properties = ±2t/Ao and = 4 /Aq, we 

find 


dUfi(E) + f2(E,v’o)) = 


8 

4 

dE{d+d-) 

\ * 

A2 ^ 


d+ 

M+P A 

8 

\d+? 

( . , id+^ 

- + - 


dE{d+d-) — 9|;[r+r_ -I- r'j^r'_ — 2d+d-\ 

O / _7 _7 \ I 2 -I 


1 dE{d+d-) 

2 d+ 


- \d+?^dE[r+r- + r+r'_ - 2d+d_] . 


We can get rid of the second order derivatives in the last term. Eirst, we note that 

dl{d-d+) = -2\dEd+\'^ + 2Red\dld+ , 
9|(r^V^^) = 2\dErl^\^ - 2-Rerl^*dlrl^ , 

and also 

Re(i!j_i9|fi+ = -\dEd+\^ - - Rer+(9|r+ 

= -lagd+p - - Rer+a|r+ , 


(S24a) 


(S24b) 


as can be found by differentiating the unitarity condition jd+p -I- = |d+p -f |r+p = 1. Combining the above two 

equations gives. 


d%{d-d+) = -4|5£;d+p - I^Br+p - ISsJ’+P - Rer;i_9|r+ - Rer^(9|r+ 
= 2|9£;?'+^P — 2Rer^^*9|;r^^ 

The second order derivatives cancel in Eq. (S23), and we are left with 


(S25a) 

(S25b) 


d%{ME)+ME,^o)) = : 


dEd++-^\ —\dEd+\ -f- 


2 , 1 |'9£;(d+d-) I 1 2 


--\d+ndEr+\^ + \dEr'+\^ + 2 \dEd+\^] 


We will next show that the r.h.s. of the above equation can be simplified, 

- \dEd+\^ + ^ + \dEr'+\2 ^ 2\dEd+\^] = -\dE{d*_r+)\^ . (S27) 

Using the unitarity of scattering matrix, = —d±r^/d'^, we have 

\dEr+\'^ = l^er+p - 4:——^lmd+dEd*_^_lmr+dEr+ + 4:y—^(lmd+dEd*_^_ f (S28) 

fo+r \d+r 

where we used ^9e{^) = in the second and third terms. Einally inserting above equation into the 

left-hand-side of Eq. (S27) and using dsid+d-) = 2zlmd^j_9£;d+, we find 


- \dEd+\^ + i + 2\dEd+\^] 

= -\dEd+\'^ E 2{lr[vd+dEd*^y - |d+p|a£;r+p -f 2Im(i+aBd!j_Imr+(9£;r+ - \d+\^\dEd+\^ 

= -\r+\^\dEd+\^ - Id+H^^r+p -fi 2 (Re(d+a£;d+)Re(r+a£;r+) - Imd+a£;d:!;_Imr:!i_(9£;r+) 

= -\r+\^\dEd+\^ - \d+\^\dEr+\’^ + 2Re{d+dEd\){r\dEr+) 

=-\dE{d*_r+)\\ (S29) 
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In the third line we used the complex number identity (Imw)^ = |u;p — (Rew)^, and also Re{d+dEd’^) 
from unitarity. 

With the above identity, we can write 


d%{h{E) + ME,^o)) 


8 

\d^ 


^ dEd+ + 


id+ 

Aq 


2 


\dE{d*_TE)\^ 


Inserting this into Eq. (S22) leads to Eq. (SIO). 


-Re(r^5£;r+) 


(S30) 


SM4. THE SCATTERING MATRICES 

The Andreev and normal scattering matrices Sa and Sn were already defined in Sec. SMS, in Eqs. (S12) and (SI5). 
In this section we express them in terms of the microscopic parameters of the model Eq. (3) in the main text. This 
is done by using the transfer matrix. The Andreev and normal scattering are discussed in Secs. SM4 A and SM4 B, 
respectively. Finally, in Sec. SM4 C we derive the normal reflection amplitudes in the Born approximation for a weakly 
reflecting junction. This result will be put to use in Sec. SM6. 

Since Eq. (3) of the main text, has only one spatial derivative, its eigenstates can be conveniently written 

with the help of the transfer matrix. An eigenstate of energy E can be written as 

(^e{x) =T{x,xq\E)<^e{xo) , (S31) 

where the transfer matrix is 

r i 

T{x,Xo] E) = T exp / dx' -TzCTz {E + — V{x')]tz — M{x')ax — ^o{x')[tx cos (p{x') — Ty sin (p{x')]) , (S32) 

Jxo 

and the ordered exponential is defined as, 

T exp f dx7(a;') = lim Aa; = Xn = xo + nAx , (S33) 

Jxo ^ 

and satisfies dxT exp dx'f{x') = f{x)T exp dx'f{x'). 

From the transfer matrices T(—oo,0), T{0,L), and T(L,oo), we can read off the scattering matrices Sa,l, Sn, 
and Sam, introduced in Sec. SM3, Eqs. (S12) and (S15). (Hereafter, we shall mostly suppress the explicit E in the 
arguments of the transfer matrix and scattering matrix.) We will start by considering the superconducting regions, 
a; < 0 and x > L. 


A. Transfer matrices outside the junction 


We assume that M{x') = 0 in the superconducting regions. The transfer matrix (S32) is then diagonal in spin- 
indices. The term iaz[n— V(x')]/v in the exponential in (S32) commutes with the other terms, and can be pulled out 
to factorize the transfer matrix. Also, in our model (3) of the main text, Aq and (p are piecewise constants so that 


T{x,xe) = e 


_ ^cr^{n-T)K,(x-xo) 


KV V 

where either a;, a:o < 0 (in which case we set (/? = 0 in h), or x, xq > L (in which case (p = tpo). Using Eq. (S34) in 
Eq. (S31) with {x, a;o) —>■ (0, —oo), and (x, Xg) —> (oo,L), leads to Eqs. (S12), (S13). 

Using the identity 


(Tz (n-T IKIX 


-Xo) _ 


1 

2 L 


(1-I-(Jz(h • r)) + {1 - az{n ■ t)) e 


(S35) 


in Eq. (S34), and requiring that $£;(±oo) = 0 we hnd the equations 

{l + az{n-T))^E{L) =0, {I - az{n- t))^e{0) = 0 , (S36) 

which lead to Eq. (S14). Combining the above three equations in Eq. (S31) shows that the wave functions decay 
exponentially outside the junction, 


~ 1 e-'"- , a; < 0 . 


(S37) 


The decay length is k ^ and approximately equal to ^ at energies E <C Aq. 



7 


B. Transfer matrices inside the junction 

In the junction we have Ao(x') = 0 in Eq. (S32). The transfer matrix T{L, 0) is therefore diagonal in the particle-hole 
components. The diagonal elements are explicitly 

T± =T exp / dx'—-(Jz (E ± [/i — V{x')\ — M{x')ax) ■ (S38) 

Jo ^ 

Here corresponds to the electron/hole transfer matrix. For brevity, we suppress their position arguments in this 
section. To relate the components of the scattering matrix to those of the transfer matrix, we write Eq. (S31) across 
the junction: 


/ 

uiW 

vi{L) 

V V^{L) 


T+ 0 
0 T_ 


wt(0) \ 

wt(0) 

n(o) 

wt(0) / 


T± 


T±i^ ) 


Comparing Eq. (S39) to the definition of the junction scattering matrix Sn, Eq. (S15), we find 


( d'+ 

[d+ 4 


1 

T+4.4 


T 




1 

T+U 


-T. 


+n 



(S39) 


(S40) 


In deriving these expressions, we used the pseudounitarity, Tj^cr^Tj- = <Tz, of the electron and hole transfer matrices. 
One can check that the pseudounitarity implies unitarity, S'±5'± = 1, of the scattering matrices. 

Particle-hole transformation P = UyXyK [K denotes complex conjugation) relates the electron transfer matrix 
r+(E) to the hole transfer matrix T-{—E) of the opposite energy. To find how the transfer matrix transforms 
under P, we can act on Eq. (S39) with P from the left. Comparing to the same equation (S39) written for the 
opposite-energy, we find 


ayTl{E)ay = T_{-E). 


For the scattering matrix, Eq. (S40), we find the relation (in both abstract and in component form) 


(S41) 


-azSliE)az = S.i-E), /J^\e) ={-E), d‘i\E) = d^^^*{-E). (S42) 

The transfer and scattering matrices are also constrained by the symmetry of the junction Hamiltonian (3) (main 
text) under the combined reversal of time and magnetic held. The time-reversal operator is T = —iayK, while Rm = 
(Jz has the sole effect of reversing the Zeeman term, as can be checked from Eq. (3): RMM{x)axRJ^ = —M{x)(Jx. 
The combined symmetry operator TRm = cf^K commutes with the junction Hamiltonian. The effect of the combined 
symmetry on Eq. (S38) is 


r± [M{x)] = axTl [M{x)]ax • (S43) 

[The square brackets remind us that T± is a functional of M{x).] By using Eq. (S43) in Eq. (S40), we hnd that the 
left and right transmission amplitudes are equal, 


d± = . (S44) 

Note also that the sole action of time-reversal on Eq. (S38) yields 

T± [M{x)] = ayTl [-M{x)]ay . (S45) 

Combining Eqs. (S43) and (S45), we see that the diagonal elements of the transfer matrices are even functionals of 
M{x), while the off-diagonal elements are odd. For the scattering matrices, Eq. (S40), it implies that the transmission 
amplitudes are even in M{x), while the reflection amplitudes are odd. 



C. Perturbative treatment of M{x) 


Because the helical edge states of a TI are protected by time-reversal symmetry from elastic backscattering, the 
junction is reflectionless in absence of M{x). In this section we consider a weakly reflective S-TI-S junction and derive 
the scattering matrix to first order in M{x). As explained at the end of the previous section, we expect the reflection 
amplitude to vanish linearly, ~ M, for small M, while the transmission amplitude remains non-zero, = 1. 
The corrections to the latter will be of order and we neglect them in this section. 

To first order in M, we can expand the transfer matrix as T{xj,Xj-i) = T^^'>{xj,Xj-i)+6T{xj,Xj-i). Here Xj, Xj-i 
are any two points in the junction. The zeroth-order-in-M transfer matrix is diagonal in spin-indices. From 
Eq. (S32) we find 

Xj_i) = exp - / [Etz^ — V{x')\azdx'. (S46) 


The first order correction 6T(xj,Xj-i) ~ M can be written in a simple form for close enough points Xj, Xj-i. 
Assuming M{x) does not vary significantly over a short interval [xj, Xj-i], we find from Eq. (S38) to lowest order in 

{Xj - Xj-i) 


ST{xj,Xj-i) = -Tz(JyM{xj){xj — Xj-i). 


(S47) 


Consider next the full transfer matrix T{L,0) of the junction. It can be decomposed as a product 


r(L,0)= lim T{xn,xn-i) ■ ■ ■T{xj,Xj-i).. .T{xi,xo), Xj = jL/N . 


(S48) 


Expanding Eq. (S48) to first order in 5T yields 

T{L, 0) = (L, 0) -f (L, 0), 

where 

N 

T*-^^(L,0)= lim T^'^\xN,Xj)ST{xj,Xj-i)T^^\xj-i,Xo) 


j=i 


- f dxT^^\L,x)Tz<JyM{x)T^'^\x,Q). 
X Jo 


Using Eqs. (S40), (S46), (S49), and (S50), we find the scattering amplitudes to lowest order in M, 


and 


i 

d± = d'. = exp [E ± (g, — V(x'))]dx ', 

V Jo 


1 ‘2i'L 

r± =-i- / dxM{x)exp — / [E ± {n — V{x'))]dx', 


(S49) 


(S50) 


(S51) 


1 Qii 

r'j. =-i- / dxM{x)exp — / [E ± {g, — V{x'))]dx'. 

V Jo X Jz; 

Alternatively, the right reflection amplitude can be written as 

1 2i f'^ 

r'^ = —i- / dxM{L — x)exp — / [E ± — V{L — x'))]dx'. 

X Jo X Jq 

This shows that r = r' if both M and V are reflection symmetric, M{x) = M{L — x) and V{x) = V{L — x). 
In a simple example of M{x) = MQ{x)Q{L — x), V = 0, and M v/L, we find 


(S52) 


(S53) 


r± = r'+ = M 


1 _ gf 

2[E±gL] 


(S54) 


as quoted in the main text below Eq. (10). 

The equations (S51)-(S52) are the main results of this section, and will be put into use in Sec. SM6. 
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SM5. SYMMETRIES OF THE HAMILTONIAN (3) OF THE MAIN TEXT 

In this section we discuss the symmetries of the model (3) of the main text which lead to degeneracies at special 
values of the phase difference (pQ. We will not consider zero energy excitations in this section - their existence was 
discussed below Eq. (S9). The Hamiltonian (3) antzcommutes with P = (TyTyK {K denotes complex conjugation). 
This particle-hole symmetry does not lead to any degeneracies in the excitation spectrum, and we shall not discuss it 
any further here. 

We start by considering a reflectionless junction, M{x) = 0 [see remark below Eq. (Sll)]. At M{x) = 0 the system 
is time-reversal symmetric (TRS) for real values of the order parameter i.e., at i^o = 0) At these values 

of the phase, the time-reversal symmetry operator, T = layK, commutes with the Hamiltonian (3) of the main text. 
Since = — I, this leads to Kramers degeneracies at (/Sq = 0 or tt, regardless of p or V{x) (see Sec. SM6 A). 

A non-zero M(x) breaks time-reversal symmetry and for generic p and V(x) this leads to the opening of gaps at 
(po = 0 and tt. Curiously, the degeneracies at po = tt are preserved it p = V(x) = 0 and M{x) is reflection symmetric 
about the junction center, M{x) = M(L — x). The degeneracy exists because there are two operators, = (TxTx and 
S 2 = TyR, that commute with at po = Here R is reflection operator, Rf{x) = f{L — x). Since these two 
operators Si and S 2 do not commute there is a degeneracy at po = tt. This degeneracy can be lifted by ^ 0, see 
Eq. (11) of the main text and the associated footnote. 


SM6. WAVE FUNCTIONS AND SPECTRUM IN A LONG ALMOST BALLISTIC JUNCTION, L > ^ 

In this section we derive the sub-gap energy spectrum and the transition matrix elements in a weakly reflecting 
(|r£^| <C 1), long {L C) junction. We shall focus on (/?o ~ tt where, as discussed in the main text, the energy 
spectrum of the high Andreev levels (E„ > Em) is characterized by avoided crossings. 

We start, in Sec. SM6 A, by deriving the spectrum and wave functions for a reflectionless time-reversal symmetric 
junction (allowing for scalar disorder). In Secs. SM6B and SM6C we derive the lowest order (oc r) backscattering 
corrections near p^ = tt to wave functions and spectrum of states n > 1 and n = 0 (Majorana, Ai), respectively. In 
Sec. SM6D we show how to use the derived wave functions to find the squared transition matrix elements, Eq. (12) 
of the main text. Finally, in Sec. SM6E, we calculate the correction oc to the slope dEM/dpo of the lowest level. 


A. Spectrum and wave functions in a reflectionless junction, r± = r± = 0 

For a reflectionless junction, Eq. (S16) factors into 

_ d+{E)d-{E)^ _ d+{E)d-{E)^ =0, — , (S55) 

which has degenerate solutions at pn = Q and tt (as expected from TRS, see Sec. SMS). Using Eq. (S51) we have 
d+d- = exp2iAL/v and Eq. (S55) takes the form [S6] 


EL jv — arccos = ± ^ -I- 7rm, to G Z . 

Aq 2 

Assuming E <C Aq, we can expand arccos A/Aq « 7r/2 — E/Aq and find a solution 

l + (-l) 




2(i + 0 


TT n + 1 — 


+ (-ir 


2(A + 0 


\Ph-tt\, neZ. 


(S56) 


(S57) 


Our assumption is valid when L ^ = v/Aq, and n L/^. The superscript in E^^^ indicates zeroth order in 
reflection amplitudes. The factor 1 -1- (—1)" vanishes for odd n and is there to make sure that levels n and n + 1 
(with odd n) are degenerate at po = tt. The states with odd (even) n have positive (negative) slope dEn ^/dpo for 
Po < TT. In the reflectionless junction we consider here, the slope changes sign discontinuously at po = tt. Disorder in 
the junction makes this switching continuous as we see below in Sec. SM6B. 

The wave functions for Andreev levels n > 1 in the reflectionless limit can be found from Eq. (SI7e), by employing 
Eq. (S16) when carefully taking the limit r —)■ 0. Consider first po < tt. For even n we have d+d_ = and 
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therefore = 0 directly from Eq. (S17e). For the odd n we find = 0. For any we find at x = L, 


0 

IveniL) = 0(7r - ipo)Vn^iL) \ ® ^ 

1 


/ 

V 




f 


J 


$ 


(0) 


odd(-^) = - Vo)UntiL) 


1 

0 

0 


+ 0((/?o - 7r)u„^(L) 


+ 0 (v?o - 7 r)v„^(L) 


1 
0 

g-ia(-E^“^)g-*Vo 




n even, (S58a) 


0 

^iipo 

0 

1 


odd. 


(S58b) 


The wave functions at any position can be constructed from the boundary values by using the 

transfer matrix introduced in Sec. SM4. In the proximitized regions they decay exponentially according to Eq. (S37). 
Inside the junction the wave functions are extended, 


/ 


0 


Ln(a;) = 0(71- - Vo)Vn^{L) 


- TT)Un^{L) 




V 


V 


) gi(^o g- I [E^'''+^^-V{x')]dx 

0 

0 

0 


0 < X < L , 


(S59) 


here given for even n with odd n has similar form). Since the Hamiltonian Eq. (3) in the main text, 

commutes with Sz, the states $„ are also S'z-eigenstates. Normalizing the wave function, J (x) = 1, 

leads to 


UntiL) = Vn^{L) = 


-\J 2{L + K„ 


K„. = -\/Al- E^. 


(S60) 


When E <C Ao, we have k ^ and = u„^(L) = Ijyj2{L + ^). 

Consider next the lowest Andreev level, n = 0. Assuming E <C Aq, we can use Eq. (S57) with n = 0 and take the 
positive solution = E^^ > 0. We find 




The wave function of the state E^ > 0 is discontinuous at ipo = tt: 


<iL) = 


\/2{L + 


/ 


0(7r - cpo) 


0 

0 

1 




- 0 (v?o - tt) 


1 

0 

^-ia{E^J^)^-i(po 


(S61) 


(S62) 


Note that in Eq. (S62) the wave function at i^o < tt is related to that at i^o > tt by a particle-hole transformation 
P = (JyTyK and simultaneous E^ —)• —E^. 


B. Backscattering corrections to spectrum and wave functions of Andreev levels n > 1 

Let us next consider the correction to spectrum near ipo = tt due to a small reflection amplitude r. The spectrum 
equation (S16) can be written as 


1 - 


d+d-e 


— 2icx 




M-l 


1 - 


d+d-e 


— 2ict 




M+P 


— e 


— 2ioc 


{d+d-r'* - r'+E^°) {d+d-r'* - r'_e-^‘^°) = 0. (S63) 



11 


(Here the scattering amplitudes and a are of course energy-dependent; we shall mostly omit the argument E from 
these quantities in this section.) 

In deriving Eq. (S63) we used the identities (derived from unitarity of the scattering matrix) r± = —d±r'^/d'^ 
and d±d'j. — r±r'^ = d±/d'^, and d± = d'^, Eq. (S44). In the last factors in Eq. (S63) we used |(i±| = 1, neglecting 
terms of order r^. We shall perform a double expansion of Eq. (S63) in both r and cpo — ir to find the backscattering 
corrections to spectrum near the degeneracy point of the clean limit. In the first two factors in Eq. (S63) we can 
replace ~ I since the corrections oc to the norms do not contribute at zeroth order in tpQ — tt. In the last two 

factors which are both oc r', we use the zeroth order vo — t: expressions, = —1 and d+d- = —. 

The former is obtained from Eq. (S55) and is the energy of the reflectionless junction at (/?o = tt, see Eq. (S57). 

Upon doing these approximations and expanding the first term in — tt, Eq. (S63) becomes 

(1 -f + [pQ-irf + \r'_ - e^*“r'+=0. (S64) 

Let us denote 6E ~ max(r, \pq — 7r|) the small correction to the energy, E = E^^'>{tt) -|- 5E. Erom Eq. (SSI) we have 
c?_i_c?_ = negligible corrections oc to the norms |(i±|. Using a{E) « a(E^^^Tr)) — SE/Aq 

and the zeroth order the solution g 2 iE(°'>(Tr)L/v _ Expanding 

the latter to first order in 6E we can solve for it in Eq. (S64): 


SE = ± 


2(i + ?) 


y {po - 7r)2 + |U - e2-r;* . 


(S65) 


The spectrum near the avoided crossing is then (see Fig. 2a in the main text) 

Eni^o) = ^ ~ 2 ) ~^ *'~^^" 2(L + g) ^^‘^° ~ (S66) 

and reproduces Eq. (S57) in the reflectionless limit. Equation (S66) is Eq. (10) in the main text. The width of the 
avoided crossing region is. 


Spn 


— e 


I (,r) 


(S67) 


The above equation is written in terms of the right reflection amplitudes r(|_. By using the identity r'^ = —d±r^/d'^, 
and d+d- = —at energy E^^'>{Tr) one can show that Eq. (S67) can be also written in terms of the left re¬ 
flection matrices, 5pn = k-i-~ ^o lowest order in U„/Ao, we have = —1 and we get 

Spn = \r+ + 1 (ti introduced in the main text. In a simple example of M(x) = MQ{x)Q{L — x),V = 0, and 

M, /X <C u/L, we find from Eq. (S54) 


(r+ -I- T_) 


v/Le^u\tt) 


(S68) 


Using the above equation in the expression for 5pn in the limit <C Aq leads to Eq. (11) of the main text. 
Next, we derive the wave functions <i>„(L) near the avoided crossing. From Eqs. (S17e)-(S17g) we find 


^n{L) 


Un^{L) \ 

e-*“e-*‘^'>M„^(L) I ’ 

Vn1;{E) / 


Unt{E) 


-d+r'l/d*_ 

e*“ - d+(l/dl)e-*“e-*‘^o ' 


As before, we expand Un^{L) in small po — t: and r'^ by inserting d+d-e We find 


(S69) 


Un^{L) 




{r'l - e-2“rV)^^^(o) 


(tt) 


{po - tt) - (-1)” V(V5o - 7r)2 -h {5pn)‘^ 




(S70) 


up to corrections of order {6E)'^. The correction to <I’Jj(L)$„(L) by r is of second order, and beyond our accuracy. 
Therefore, the normalization factor u„^(L) is determined by requiring $jj(L)$„(T) = 1/(L-|-^), see Sec. SM6A. This 
yields 


- 1 1 _ 1 1 ( icn (7^0-^) 

\/2{L + ^) Vs y \/{‘fo ~ + (^T^n)^ 


Writ (A) 


(S71) 
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Inserting (S71) into Eq. (S70), we find 




^/2(E + 0^/2 


11 + (-ly 


(<Po - tt) 


- 7r)2 + {SlfinV ’ 

^0), 


(S72) 


where 'dn = arg(r^ — e and all the energies are evaluated at = En^ir). It is convenient to introduce 


so that 


smt/n = 

‘^n oddt(-^) “ 
event ~ 


a + 


{ipo - tt) 


t2 Y Vivo - 


COSf/n = 


d- 


((/?o - tt) 


V^Y ^((^0 - 7’')2 + (( 5 (/ 9„)2 ’ 


sin 

■\/2(L + C) 

cos 6i„ 

V‘^{L + C) 


«noddT(^j- y2(LTe) 

Zi„evenT(7.)-Ze 6 ^^(L + e)' 


(573) 

(574) 

(575) 


Finally, in terms of the wave functions (S58a), (S58b) of the reflectionless junction, the perturbed wave functions (S69) 
near the avoided crossing are 


$(o) (L) = 

n evenx / 


cost 


VWTl) 


1 




( 1 

\ 

0 



sin On 

V2(A + 0 

0 

^-ia{En) 

g-*Vo 

V 1 

) 


\ 0 

/ 




VW+l) 


/ 1 



( ® 

\ 

0 


sin On 

^-ia(En) pivo 


^-ia{E„) 

p-iVo 


0 


V 0 

) 


1 1 

/ 


even, (S76a) 


n odd. (S76b) 


C. Backscattering correction to the wave function 4 >ai of Andreev level n = 0 


In this section we will calculate the backscattering correction 5^m to the wave function of the Majorana doublet 
(n = 0 state). The presence of a Ej^ = 0 state at = zi" is not removed by any r. An illustration of this is that 
5ipn=o = 0 in Eq. (S67), as can be verified by using the properties r)|_(0) = —r)j‘(0) [see Eq. (S42)] and = —1 

in that equation. Unlike the crossing of the levels {E_\ 4 ,—Em), the degeneracy at (^o = 0 of the pair {Em,Ei) is 
lifted by r. A finite value of r modifies also the slope dEj^/d(po at (/?o = tt by a small correction oc r^, calculated in 
Sec. SM6E. Since we are interested here only in the first-order corrections, we will neglect the modihcation of E_v(, 
and focus on the first-order correction to the wave function. 

Let us first consider phase (/?o < tt. Then the unperturbed wave function (S62) has V^^{L) = 0. From the second 
equation in Eq. (S69), we have to lowest order in 


(A) 




(S77) 


Both the numerator and denominator in the above equation vanish at (po = to as can be checked by using the properties 
at i?_v( = 0: r[|_(0) = —r[j'(0), = i, and d+d- = 1. Taking the limit (po ^ tt yields (we use dfdipo = d^^EdE+d^g) 


UMtV) 


, 2d^gEM{-r+dVdEd+ + -b 

* -2d^gEMidVdEd+ - idEaiE))+i 


(S78) 


To zeroth order in r the transmission amplitude satisfies dVdsd^ = iL/v, see Eq. (S51). Also, for ipo 7r“ we have 
—2d^pgEM = ot/{L + y). Finally, using dEC({E) « —and the condition r(|_ = —following from unitarity of 
S+, we get 




V 

2L 


d\dEr\ 


1 

VVE+V) 


Vo < TO, 


(S79) 
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where we also used v_\ 4 -f(L) = Xj + ^). The corrected wave function at (/?o < is 


where 




1 

V2(i + e) 



( ^ 


( 1 ^ 





0 ' 



0 1 

+ S<t>M 




1 1 } 


\ 0 J 



V30 < 71", 


= ^d\dErX , 


(S80) 


(S81) 


and the scattering amplitudes and their derivatives are evaluated at ill = 0. As seen in Eq. (12) of the main text, the 
correction modihes the transition matrix elements slightly. 

The wave function at i^o > tt can be found by solving for vq^{L) in Eq. (S17e) with the help of (S16), and then 
doing a calculation similar to the one above. A more convenient way is to use the particle-hole transformation P [see 
remark below Eq. (S62)] on Eq. (S80). We find 




— , 8(7r — (Pq) 

y2(LTO 

— , Qj^Po — Tt) 

y2(LT0 



/ 0 \ 

/ 

1 ^ 



^-ia{EM) giVo 


0 



0 


^-ia{EM)g-iVo 



V 1 


0 



f 1 

\ 

f \ 



0 


^-ia{EM) giVo 



^-ia{EM)g-iVo 


0 



V 0 


V 1 / 



(S82) 


In the simple example of M{x) = MQ{x)<d{L — x), V = 0, and M <C v/L, we find from Eqs. (S51), (S54), (S81) 


S<X>M 


ML /I 
~ 1^2 



(S83) 


This result was quoted below Eq. (12) in the main text. 

The example used in the main text in the context of a NW was a (5-impurity at the NW-S interface, M{x) = 
V(x) = UoS(x — L). The scattering matrices for a more generic 5-impurity at x = a;o are calculated in Sec. SM8. In 
our specihc example we hnd to lowest order in uq/v, 


= 


Uq 

V 


( 



V 


(S84) 


D. The transition matrix elements 


In the main text we saw that the oscillator strengths appearing in the admittance Y are given by the matrix 
elements of the perturbation From Eqs. (3) and (6) of the main text we can read off 9^g'H*'°^(x) = 

—0(x — L)Ao[ra; sin po + Ty cos po] and 


= -Ao 


dx^'j^{x)[Tx sin(,i)o + Ty cos (,i)o]d>n(x) . 


(S85) 


The integration domain covers only the proximitized region, where the wave functions decay exponentially from their 
boundary values at x = L, see Eq. (S37). Using the wave functions at x = L near po = tt [derived above, Eqs. (S76a), 
(S76b) and (S82)] in Eq. (S85), we hnd the transition matrix elements to lowest order in reflection amplitudes and 
\po — 7r|. Squaring the matrix elements leads to Eq. (12) of the main text. 


E. Backscattering correction to the slope dEj^/dpo of Andreev level n = 0 


In this section we calculate the terms denoted (!l(|r+|^) in Eq. (SIl). To access these second order corrections, we 
need the transmission amplitudes d± to that order. They can be obtained by expanding the transfer matrix to second 
order in M. One hnds 


d±{E) = 4(A) = (1 - i|r±(A)|2)4°)(A), 


(S86) 
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where r±{E) is given by Eq. (S52) and = exp ^ J^[E ± {p — V{x'))]dx' from Eq. (S51). The form (S86) 

satisfies the requirement |d±P + |r±P = 1 to second order in r. 

We can now insert Eq. (S86) in the denominator in Eq. (SIO) and keep terms up to second order in r. By using 
the properties dEd^p = Slid 


|5£;d++ifi+Ao^ - dE{d+r*_) 


L+l 

V 


(l-|r+n- 


iL 

OEr+ -r+ 

V 


&% E = 0, and 


d^poEMM « ± 


1 V 


1 , 


1 dEr+ - 


2 L + M ^ 2 '^+''+2 (L + 0 " 


E=0 


(S87) 


(S88) 


This is the equation referred to in the footnote 30 of the main text. In the simple example with M{x) = M<d{x)Q{L — 
x), E = 0, and M,^ ^ v/L, we have dEr+ = to lowest order in plipjE). The second term in (S88) is then 
negligible, and we find 


dp,Ej^{E) « ±2;^(1 - ^k+(0)n = ±^^M+(0)l • 


(S89) 


Note that |r+(0)| = |r_(0)|. In the second equality we used Eq. (S86). 


SM7. EXTENSION TO NANOWIRE JUNCTION 


The effective Hamiltonian = f cia;4tf'H*'°^'I'/2 of the nanowire junction is [S7, S8] 

{x) = -Tz ^ - iaTzdzdx - pTz + U{x)tz + Bdx + A(x) [tx cos p{x) - Ty sin p{x)] , (S90) 

A(a;) = A[0(—x) + 0(a; — L)], ip{x) = poQ{x — L). (S91) 

Here m is the effective mass of the electrons, a is the SO coupling constant, U{x) is a scalar disorder potential, B is 
the Zeeman energy, and A is the proximity-induced gap. The Pauli matrices d act on the spin-space and the Nambu 
spinor is “V't)- Coupling to a weak microwave source is treated as in the main text. 

We will next show that at low energies the spin-orbit coupled nanowire is very similar to the TI edge, and its 
low-energy Hamiltonian is identical to Eq. (3) of the main text. Eor illustration, we shall focus on the limit of large 
Zeeman energy, B ^ ma^,E, and set /x = U{x) = 0. Near the Fermi points k « ±kz = E\/2mB, the effective 
Hamiltonian is El^ = f dx(j)^li}^g(j)/2 with [S9] 

"Hgg {x) = -ivTzUzdx + Ao(x)[tx cos ip{x) - Ty sin (p{x )], (S92) 

where v = kzlm, Ao(a;) = A{x)2ma/kz, and = {R\ L\ L, —R) is the slow mode field operator. The Pauli 
matrices <Jx,y,z and Tx^y^z act on the respective spaces of left and right movers, and particles and holes. The above 
Hamiltonian is derived from Eq. (S90) by using the projections [S9] 

^^{x) « - 7+e-*"""T(x), 


and assuming slowly varying operators L{x) and R{x). Here 7± = (1 ± malkz) j\/2. 

The low-energy Hamiltonian (S92) has the same form as the Hamiltonian (3) of the main text with p = V = M = 0, 
and therefore the bound state spectrum of the clean {U = 0) nanowire junction is the same as that of a reflectionless 
S-TI-S system (with of course different Dirac velocity v and pairing gap Aq). 

Above we derived Eq. (S92) for a clean nanowire with no potential disorder U{x). For a simple model of disorder, 
consider a d-impurity at position a; = xq in the junction, 

Himp = X] / dxuQ5{x - xojx/'aV’a . (S94) 
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or U{x) = Uq5{x — Xq) in Eq. (S90). Using Eq. (S93), the projection of -ffimp to low energies is 

-ffimp = Mo[7- + 7+] + -^^'''(a^o)^(a:o)] (S95) 

+ mo27+7_ [e-2*'=^"»i?^(a:o)U(a;o) + e2*'=^"»Ut(a:o)i?(a:o)] • 

The phases in the second term can be gauged away by setting R{x) —>■ e~‘^^’^^^°R{x). The projected i?imp leads to an 
additional term of the form 


in Eq. (S92). Here 


'^imp = + moax]6{x - xq) , 


Vo = Mob- + 7+], 
Too = 2^07+7- , 


(S96) 


(S97) 


are the strengths of forward and backscattering due to i7imp- The effect of on the sub-gap spectrum is studied 
in Sec. SMS below. Generally, the presence of an impurity leads to lifting of degeneracies at ipo = 0, n (without 
removing Ejotln) = 0 of course). In particular, we find that an impurity at the NW-S interface, Xg = L, results in a 
width 5ipn oc rag/v of the avoided crossing at (pg = tt. This result was quoted in the main text, Eq. (13), and will be 
derived in the next section. 


SMS. SINGLE 5-IMPURITY - SCATTERING AMPLITUDES AND EQUATION FOR SPECTRUM 


For a single 5-impurity the transmission and reflection amplitudes can be calculated explicitly for any impurity 
strength. Consider an impurity 


M(x) = mg5{x — Xg), V{x) = VgS(x — Xg) . 


(S98) 


where 0 < Xg < L. We have encountered such an impurity at xg = L in the main text in the paragraph of Eq. (13) 
(see also Sec. SM7). In this section we will study the effect of the impurity on the spectrum of Andreev levels and 
derive Eq. (13). 

For the above 5-impurity the junction transfer matrix can be written as 


T(L, 0) = T(L, 0). 

where T{L,xo) and r(xo,0) correspond to reflectionless propagation, 

T{L,xo) = eH^-^o)[Er^+>^W. ^ T(xo,0) = 

Inserting Eq. (SlOO) into Eq. (S99), we find in the space of left and right movers 

-lL[Er,+^]l 


jT e“«— e 


r(L,o) = 

The complex numbers z, w satisfy |wp -I- jzp = 1. Explicitly they are 


1 


Mo 


I COS-Rvl-mlEi 

V '' - m, 


1 


-1 


2 sinyV^o 
0 


w = z 


mo 




1 

sm -\/Vn — m. 


0 ■ 


Using Eq. (S40) and above T{L,0), the normal state scattering matrices are 

^lL[E+ii\^ 


r+ 



d+ 

r'J-[ 


r_ 

d'_\ j 


d- 

r'_ ) = i 



W 




(S99) 

(SlOO) 

(SlOl) 

(S102a) 

(S102b) 

(S103a) 


(S103b) 



16 


We see that |r£^| 


|w| and |ci±| = \z\ are independent of energy. The above-derived transmission amplitudes satisfy 


^ ^ _ ^2LE 

d*_~ d%~ 


(S104) 


With the above properties, the equation for spectrum, Eq. (S16), simplifies significantly. We find [SIO] 


I 7i 2 / 7-i\i I i2 ^{L — 2xo)-E' „ 

|dp cos ifo — cos 21- oi{E)] — \r\^ cos-= 0, 

V 


where we denote |r£^| = |r| and |(i±| = |(i|, and 

\d? = 


1 




= 1 - Ir| 


(S105) 


(S106) 


If the impurity is at the interface, xq = 0, L, a gap of order \r\{v/L){^/L) opens for the high Andreev levels at 
(po = TT (for simplicity, we assume |r| <C 1 here). Equation (S67) gives the width of the resulting avoided crossing 
region. 


Sipn 



v/L 


(S107) 


If the impurity originates from a scalar impurity in a NW junction (see Sec. SM7), we can relate |r| to the microscopic 
impurity strength uq. This is done by using Eq. (S97) to express vq and mo in terms of the microscopic parameters, and 
using them in Eq. (S106). We find that when ug <C V\/B jmd^^ the reflection coefficient becomes |rp = (4 -|- 

with Z = 2mofv « 2uq/v. 

For weak backscattering (ug v) we insert |r| « |uo|/u <C 1 in the expression for Sipn to find Eq. (13) of the main 
text. The phase in this case is = {r'* — = (uo/|'Uo|)(l — iug/v). 

Finally, we note that an impurity in the middle of the junction does not lift any of the degeneracies at (fg = tt. 
This can be seen from Eq. (S105) where setting xg = L/2 yields 


cos 2[-«(-£')] = MP cos (fig — |rp . (S108) 

Setting ifg = TT and using |rp = 1 — |(ip makes the right-hand-side independent of |r|. The spectrum a,t (pg = tt 
is therefore identical to that of a reflectionless junction and has degeneracies. These degeneracies arise because of 
destructive interference between the scattered electron and hole optical paths when the impurity is exactly in the 
middle of the junction. This result is in agreement with Ref. [SIO]. 
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